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We study in this paper the resistivity encountered by Ising itinerant spins traveling 
in the so-called Ji — J2 frustrated simple cubic Ising lattice. For the lattice, we take 
into account the interactions between nearest-neighbors and next-nearest-neighbors, Ji 
and J2 respectively. Itinerant spins interact with lattice spins via a distance-dependent 
interaction. We also take into account an interaction between itinerant spins. The lattice 
is frustrated in a range of J2 in which we show that it undergoes a very strong first-order 
transition. Using Monte Carlo simulation, we calculate the resistivity p of the itinerant 
spins and show that the first-order transition of the lattice causes a discontinuity of p. 
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1. Introduction 

One of the most fascinating subjects in condensed matter physics is the study of the 
resistivity encountered by conduction electrons in crystals. Fifty years ago, the effect 
of the magnetic ordering of the crystal on the electron resistivity began to attract 
investigations. It has been shown that at very low temperature (T) the resistivity 
is dominated by spin-wave scattering, the spin resistivity p is then proportional to 
in ferromagnets."'^ ^ In the region of the ferromagnetic phase transition, p shows 
a peak at the transition temperature Tc similar to the magnetic susceptibility. De 
Gennes and Friedel^ have suggested that this behavior is due to the spin-spin cor- 
relation. Several appro xima tions have been used to treat this correlat ion a ppearing 
in several formulations jn particular in the Boltzmann's equation.'^'^ Recently, 
we have introduced a Monte Carlo (MC) simulation technique to deal with the spin 
resistivity. Our results for ferromagnets are in agreement with other theories, in 
particular the existence of the peak at Tc and its dependence on the strength of 
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magnetic field and density of itinerant spins!^^ In unfrustrated antiferromagnets, 
there has been a few investigations. Some theories predicted the absence of a sharp 
peak at Tc-^ We have shown by MC simulations that this is true, however the 
form of the rounded peak depends on the crystal structure and other interaction 
parameters."'^'-' -'^ - 

Experimentally, there has been a large nurnber of works dealing with th e spin 
resistivity in different magnetic compounds.t^ ^^ l ^^ l ^^ l ^^ l ^^ l ^ '' l ^^ l ^^l^^ l^^ These 
experiments show the existence of an anomaly of p at the magnetic phase transition. 
The shape of this anomaly depends on the material. 

One class of interesting magnetic systems is called "frustrated systems" intro- 
duced in the 1970s in the context of spin glasses. These frustrated systems are very 
unstable due to the competition between different kinds of interaction. However they 
are periodically defined (no disorder) and therefore sub ject to exact treatments. 
This is the case of several models in two dimensioni^SI^ but in three dimensions 
frustrated systems are far from being understood even on basic properties such as 
the order of the phase transition (first or second order, values of critical exponents, 
...). Let us recall the definition of a frustrated system. When a spin cannot fully 
satisfy energetically all the interactions with its neighbors, it is "frustrated". This 
occurs when the interactions are in competition with each other or when the lattice 
geometry does not allow to satisfy all interaction bonds simultaneously. A well- 
known example is the stacked triangular lattice with antiferromagnetic interaction 
between nearest-neighbors. The frustration in spin systems causes many unusual 
properties such as large ground state (GS) degeneracy, successive phase transitions 
with complicated nature, partially disordered phases, reentrance and disorder lines. 
Frustrated systems still constitute at present a challenge for investigation methods. 
For recent reviews, the reader is referred to Ref. ■ . 

Motivated by their exotic behaviors, we have studied some frustrated systems 
and found that p depends drastically on the range of interaction, and that p shows 
a discontinuity at Tc reflecting the first-order character of the phase transition.!^ 

This work aims at confirming the fact that in systems with first-order transi- 
tions, p should have a discontinuity at Tc- For that purpose, we consider in this 
paper so-called Ji — J2 simple cubic lattice with Ising spins. This system is known 
to undergo a very strong first-order transition in the Heisenberg case.^^SI ^j-^g Ising 
case studied here shows also a very strong first-order transition as shown below. 

In section 2, we present our model and MC method. The results are shown in 
section 3. Concluding remarks are given in section 4. 



I 



I 



21, 2010 1:30 WSPC/INSTRUCTION FILE Jl-J2-niplb 



Spin Resistivity in the Frustrated J\ — J2 Model 3 

2. Model and Method 
2.1. Model 

We consider the simple cubic lattice shown in Fig. [1] The spins are the classical 
Ising model of magnitude S — 1. The Hamiltonian is given by 

T~L — — Ji Si.Sj — J2 Si.Sm (1) 

(i,m) 

where is the Ising spin at the lattice site i, is made over the NN spin pairs 

Sj and Sj with interaction Ji, while is performed over the NNN pairs with 

interaction J2. We are interested in the frustrated regime. Therefore, hereafter we 
suppose that Ji — ~J (J > 0, antiferromagnetic interaction, and J2 — —rjj where 
77 is a positive parameter. The ground state (GS) of this system is easy to obtain 
either by minimizing the energy, or by comparing the energies of differe nt s pin 
configurations, or just a numerical minimizing by a steepest descent method.123 We 
obtain the antiferromagnetic (AF) configuration shown by the upper figure of Fig. [2] 
for 1 J2I < 0.25 1 Ji I, or the configuration shown in the lower figure for | J2I > 0.25 |Ji|. 
Note that this latter configuration is 3-fold degenerate by putting the parallel NN 
spins on x, y or z axis. With the permutation of black and white spins, the total 
degeneracy is thus 6. 

The phase transition of this model in the frustrated region (IJ2I > 0.25|Ji|) 
has been studied by Pinettcs and DiejPSl in the case of the Heisenberg model. It 
has been found that the transition is strongly of first order. The ordered phase 
is very unstable due to its degeneracy. As will be shown below, the case of the 
Ising case shows an even stronger first-order transition. It is therefore interesting 
to investigate the resistivity of itinerant spins traveling across of such a system. 




Fig. 1. Simple cubic lattice with nearest and next-nearest neighbor interactions, Ji and J2, indi- 
cated. 

The model we use here to study the spin transport takes into account the fol- 
lowing interactions: 

• Interaction between lattice spins given by Eq. ([1]) 
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Fig. 2. Simple cubic lattice. Up-spins: white circles, down-spins: black circles. Upper: Ground state 
when IJ2I < 0.25|Ji|, lower: Ground state when IJ2I > 0.25|Ji|. 

• Interaction between itinerant spins and lattice spins given by 

nr = -Y.hja...S, (2) 

where cTi is the Ising spin of itinerant electron and denotes the interac- 
tion that depends on the distance between an electron i and the spin Sj at 
the lattice site j. We use the following interaction expression : 

4, =/oe-"'-- with r,,-|r--r,| (3) 

where Iq and a are constants which will be chosen in section [2.21 

• Interaction between itinerant spins: In the same way, interaction between 
itinerant electrons is defined by 

'Hm = - ^ KijaiJj (4) 

hi 

if,,, = Koe-^^^^ (5) 

with (7i the spin of itinerant electron and Ki,j the interaction that depends 
on the distance between electrons i and j. The choice of the constants Kq 
and P is discussed in 12.21 

• Chemical potential term: Since the interaction between itinerant electron 
spins is attractive, we need to add a chemical potential in order to avoid a 
possible collapse of electrons into some points in the crystal and to ensure a 
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homogeneous distribution of electrons during the simulation. The chemical 
potential term is given by 

He = DVMr) (6) 

where n(r) is the concentration of itinerant spins in the sphere of D2 radius, 
centered at r. D is a constant parameter appropriately chosen. 
• Electric field term: 

-He = -eel (7) 

where e is the charge of electron, e the applied electrical field and £ the 
displacement vector of an electron. 



2.2. Choice of parameters 

Note that the effect of the crystal magnetic ordering on the resistivity is dominated 
by the first two interactions. We will show below results obtained for typical values 
of parameters. The choice of the parameters has been made after numerous test 
runs. We describe the principal requirements which guide the choice: 

i) We choose the interaction between lattice spins as unity, i. e. \ J\ = 1. 

ii) We choose interaction between an itinerant and its surrounding lattice spins 
so as its energy Ei in the low T region is the same order of magnitude with that 
between lattice spins. To simplify, we take a = 1. 

iii) Interaction between itinerant spins is chosen so that this contribution to the 
itinerant spin energy is smaller than Ei in order to highlight the effect of the lattice 
ordering on the spin current. To simplify, we take /? = 1. 

iv) The choice of D is made in such a way to avoid the formation of clusters of 
itinerant spins (coUapse) due to their attractive interaction [Eq. ([5])]. 

v) The electric field is chosen not so strong in order to avoid its dominant effect 
that would mask the effects of thermal fluctuations and of the magnetic ordering. 

vi) The density of the itinerant spins is chosen in a way that the contribution 
of interactions between themselves is neither so weak nor so strong with respect to 
E,. 

Within these requirements, a variation of each parameter does not change qual- 
itatively the results shown below. As will be seen, only the variation of Di does 
change drastically the results. That is the reason why we will study in detail the 
effect of this parameter. For larger densities of itinerant spins, the resistivity is 
larger as expected because of additional scattering process between itinerant spins. 

We fix Ji = — J = — 1 (AF interaction) for NN coupling of lattice spins as said 
above. The energy is thus measured in the unit of J. The temperature is expressed 
in the unit of J//cb. The distance (Di, -D2) is in the unit of a, the lattice constant. 
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2.3. Simulation Method 

We consider a film with a thickness of cubic ceUs in the z direction. Each of 
the xy planes contains x Ny cells. The periodic boundary conditions are used 
on the xy planes to ensure that the itinerant electrons who leave the system at the 
second end are to be reinserted at the first end. For the z direction, we use the 
mirror reflection at the two surfaces. These boundary conditions conserve thus the 
average density of itinerant electrons. Dynamics of itinerant electrons is created by 
an electric field applied along the x axis. 

Simulations are carried out in the following manner. The lattice spins are equili- 
brated at a temperature T. Itinerant spins are then injected into the system. Before 
calculating thermal averages of transport properties, we equilibrate itinerant spins 
during a large number of MC steps. The multi-step averaging procedure has been 
used to get good statisticsff^ Averaging is made between re-equilibrating periods 
of lattice and itinerant spins to explore a maximum number of microscopic spin 
configurations. 



3. Results 

We show first the result of the lattice alone, namely without itinerant spins. The 
lattice in the frustrated region, i. e. \J2/Ji\ > 0.25, shows a strong first-order 
transition as seen in Fig. [3] The sublattice magnetization and the energy per spin 
as functions of T for J2 = — 0.26| Ji| for the lattice size — Ny — 20, — 6 
show a discontinuity at the transition temperature. To check further the first-order 
nature of the transition, we have calculated the energy histogram at the transition 
temperature Tc- This is shown in Fig. |4l The double-peak structure indicates the 
coexistence of the ordered and disordered phases at Tc- The distance between two 
peaks represents the latent heat. 

Now we consider the lattice with the presence of itinerant spins. As far as the 
interaction between itinerant spins is attractive, we need a chemical potential to 
avoid the collapse of the system. The strength of the chemical potential D depends 
on Kq. We show in Fig. [5] the collapse phase diagram which allows to choose for a 
given Kq, an appropriate value of D. 

We show now the main result on the spin resistivity versus T for | J2I = 0.26| Ji| 
for several values of Di. Other parameters are the same as in Fig. [3l As said in 
section [2?2l within the physical constraints, the variation of most of the parameters 
does not change qualitatively the physical effects observed in simulations, except 
for the parameter Di. Due to the AF ordering, increasing Di means that we include 
successively neighboring down and up spins surrounding a given itinerant spin. As 
a consequence, the energy of the itinerant spin oscillates with varying Di, giving 
rise to the change of behavior of p: p can make a down fall or an upward jump at 
Tc depending on the value of Z^i as shown in Fig. [HI Note the discontinuity of p at 
Tc- This behavior has been observed in the frustrated FCC antiferromagnet.'^ 
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Fig. 3. Upper: Sublattice magnetization M versus T, Lower: Energy versus T, for | J2I = 0.26| Ji|, 
N^=Ny = 20, = 6. 
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Fig. 4. Energy histogram taken at the transition temperature Tc for J2 = — 0.26| Ji|: black circles 
are for = Ny = 20, Nz = 6,Tc = 1.320, void circles for = Ny = 30, Nz = 6,Tc = 1.320 
and black triangles for A^^ = A^j; = 20, = 10, Tc = 1.305. Other parameters are Ig = Kg = 0.5, 
Di = 0.8a, D2=a, D = l,e=l. 
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Fig. 5. Phase diagram in the plane {Ko,D). The collapse region is in black, for IJ2I = 0.26|Ji|. 
Other parameters are Di = D2 = a, Iq = 0.5, e = 1. 
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Fig. 6. Spin resistivity versus T for IJ2I = 0.26|Ji| for several values of Di: from up to down 
Di = 0.7a, 0.8a, 0.94a, a, 1.2a. Other parameters are = Ny = 20, N^, = 6, = Kq = 0.5, 
D2 = a, D = 1, e = l. 

4. Concluding Remarks 

From the results shown above for the strongly frustrated Ji — J2 model, we conclude 
that the spin resistivity reflects the nature of the first-order transition: it undergoes 
a discontinuity at the transition temperature. The fact that as T 0, p increases 
slowly stems from the freezing of itinerant spins with decreasing T. This has been 
experimentally observed in ferromagnets and antiferromagnets as seen in Fig. 11 
of the paper by Chandra et al. on CdMnTe,[I3l Fig. 2 of the paper by Du et al. 
for MnFeGej^ Fig. 6a of the paper by McGuire et al. on AF superconductors 
LaFeAsO,'^Fig. 2 of the paper by Lu et al. on AF LaCaMnO ,11^1 and Fig. 7 of the 
paper by Santos et al. on AF LaSrMnO.H^l 

We hope that these MC results will stimulate further theoretical calculations 
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and experiments. 
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